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Bond percolation on infinite heavy-tailed power-law random networks lacks a proper phase transi- 
tion; or one may say, there is a phase transition at zero percolation probability. Nevertheless, a finite 
size percolation threshold q c (N), where N is the network size, can be defined. For such heavy-tailed 
networks, one can choose a percolation probability q(N) = pq c (N) such that lim (q — q c (N)) — 0, 

N— >oo 

and yet p is arbitrarily large (such a scenario does not exist for networks with non-zero percolation 
threshold). We find that the critical behavior of random power-law networks is best described in 
terms of p as the order parameter, rather than q. This paper makes the notion of the phase transition 
of the size of the largest connected component at p = 1 precise. In particular, using a generating 
function based approach, we show that for p > 1, and the power-law exponent, 2 < r < 3, the largest 
connected component scales as ~ N /T , while for < p < 1 the scaling is ~ N t ■ here, the 
maximum degree of any node, kmax, has been assumed to scale as N 1 ^ . In general, our approach 
yields that for large N, p 1, 2 < r < 3, and kmax ~ N 1 ^ , the largest connected component scales 
as ~ / 9 1/(3 - T) iV 1 - 1/T .Thus, for any fixed but large N, we recover, and make it precise, a recent result 
that computed a scaling behavior of q 1 ^ 3 ~ T ^ for "small g". We also provide large-scale simulation 
results validating some of these scaling predictions, and discuss applications of these scaling results 
to supporting efficient unstructured queries in peer-to-peer networks. 



I. INTRODUCTION 

Percolation on random graphs with heavy tailed 
power-law degree distributions is known to possess 
certain unique properties. Unlike most cases of finite 
and infinite dimensional percolation, the percolation 
on an infinitely large, heavy tailed power-law (PL) 
random graph lacks a percolation threshold, i.e., the 
phase transition happens at zero percolation proba- 
bility. This is the source of many interesting critical 
behavior in heavy-tailed random PL networks that 
are remarkably different from those observed in net- 
works with non-zero percolation thresholds. In this 
paper, we are particularly interested in the scaling 
of the size of the largest connected component close 
to the criticality, for networks of finite but large size 
N. To that end, we first make the notion of closeness 
to criticality precise, and provide a summary of our 
results. 



A. Summary of the Results 

Consider a heavy-tailed PL network with exponent 
2 < t < 3. A finite size percolation threshold q c (N) 
can be defined for these networks that goes to zero 
as iV goes to infinity. 
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Consider percolation at a probability q(N) = 
pq c (N). After bond percolation with probability 
q(N), a total of around H(N, p)=q(N)N{k) /2 links 
remain in the network, where (fc) is the average de- 
gree. Of these H(N,p) links, a number Q(N,p) of 
them belong to the same largest connected compo- 
nent. This paper answers the following questions: 
How does Q (N, p) scale with N and pi How does 
the fraction Q{N , p) /Ti{N , p) scale with TV and pi In 
other words, what is the chance that a link left af- 
ter the percolation belongs to the giant connected 
component? Does this probability converge to some 
constant? 

The answers to these questions are known for net- 
works with finite percolation threshold q c — q c (oo) > 
( note that p is meaningful only when p < l/q c , 
otherwise one has q > 1). In particular, Q{N, p) ~ N 
for 1 < p < l/q c , and Q(N,p) = 0(k max \ogN) for 
< p < 1, where k max is the largest degree of the 
nodes in the network. When p — 1 , Q (N, p) is known 
to scale as N 2 / 3 0] for most random networks under 
certain constraints. The probability that a random 
link after percolation belongs to the largest connected 
component, G(N, p)/Tt(N, p), converges to a constant 
Ld(p) > when l/q c > p > 1, and goes to zero at least 
as fast as fc —^ ogAf when < p < 1. 

The answers to the above questions for heavy- 
tailed PL networks turn out to be distinctly dif- 
ferent from those for random networks with finite 
percolation thresholds. In order to capture these 
differences succinctly, we define a scaling exponent, 
A T (p)4lim^ oo log(0(iV,p))/log(Ar). Fig. d shows 
a schematic of some of the results that we derive, 
based on our application of the generating function 
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approach. 

For 1 < p < N 1 ^, we show that Q(N,p) ~ 
(fc)iV 1 - 1 / r (see Eqns. HHHSJl. In other words, the 
scaling exponent is A T (p) = 1 — 1/t for any p > 1, 
and thus, the size of the giant connected component 
does not scale linearly with N for constant p > 1. 

For < p < 1, the scaling exponent A T (p) = 
1 — 2/r. This last statement follows from the ob- 
servation that after percolation, the expected degree 
of a node with deg ree k max will reduce to pqckmax- 
This last statement follows from the observation that 
after percolation, a node with degree k max will have 
an average degree of pq c k max . Taking k max ~ N 1/t 
and using the fact that q c ~ TV 1-3 / 1 " (see Eqn. Hill 
gives the desired result. 

Therefore, if we consider p as the order parame- 
ter, then there is a phase transition at p — 1 for all 
random networks (irrespective of whether it is heavy- 
tailed with unbounded variance or not) except that 
the phase transition levels are different (Fig. [TJ. This 
provides a unified view of the phase transition phe- 
nomenon for largest connected components in ran- 
dom networks. 

The above result enables us to calculate the fol- 
lowing interesting quantity; G(N, p)/H(N, p) (see 
Eqn. |20I). We find that for p > 1, the ratio 
g(N,p)/H(N,p) ~ ]\r( T_2 )/ T for 2 < r < 3. In other 
words, the probability that a randomly chosen link 
after the percolation belongs to the giant connected 
component goes to zero as N^ T ^ 2 ^ T for 2 < r < 3. 
Only for r — 2, a finite fraction of all the links that 
remain after the percolation belong to the giant con- 
nected component. 



B. Relation to Previous Work 

The finite-size scaling of the cardinality of the 
largest percolation components close to criticality has 
always been a source of interest. While for finite di- 
mensional percolation problems, the issue is still sub- 
ject to a great deal of controversy 0, it has been 
successfully resolved for percolation on most random 
graphs (that are examples of infinite dimensional per- 
colation). In particular, for an Erdos-Reyni random 
graph with percolation threshold q c , the size of the 
largest connected component is known to scale as 
8(7V 2 / 3 ) when q = q c + 0(N- 2 ^ 3 ). Thus, right at 
the percolation threshold q c , the size of the largest 
connected component scales with the network size as 
~ 7V 2/3 0]. Furthermore, it can be shown that for 
any fixed large N, the size of the largest connected 
component scales linearly as \q — q c \ when q — q c is 
small. This indicates, among other things, the conti- 
nuity of the size of the largest connected component 
as a function of q at q — q c . Note that in such cases, 
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FIG. 1: The scaling exponent, A T (p), undergoes a first 
order phase transition for all random networks at p = 1. 
Recall that A T (p) = limjv-.cxi log(Q(N, p))/ log(AT), where 
Q(N, p) is the size of the largest connected component in a 
network of size N after bond percolation with probability 
pq c . The top figure shows the transition for a heavy- 
tailed power-law random network while the bottom figure 
corresponds to a random network with a bounded k max , 
and therefore with non-zero percolation threshold. 



the percolation threshold q c is finite, and independent 
of N. In 0, this result has been extended to random 
networks on any given degree distribution, subject to 
certain constraints. 

Recent results in Q have shown that critical scal- 
ing properties of random power-law networks with 
exponents r in the range [2,4], are significantly dif- 
ferent from those of most other random networks. For 
heavy-tailed PL networks (i.e., 2 < r < 3), which is 
the case of interest in this paper, results in |J] sug- 
gest that the size of the largest connected component 
scales as \q — qd 1 ^ 3- ^ for any fixed but large N (as 
opposed to a linear scaling of \q — q c \ for most other 
random networks). The approach used in this work 
is also based on the generating function approach 
adopted in this paper. Our goal, however, is to study 
the scaling laws as a function of N for a fixed q/q c - 

We first argue that the relevant order parameter in 
heavy-tailed power-law random networks is p = q/q c 
and critical percolation properties are best explained 
in terms of p. We then derive, in a unified way, the 
scaling of the size of the largest connected compo- 
nent as a function of p and N. We show that for 
any fixed p > 1, the size of the largest connected 
component scales as TV 1-1 / 7 ", in contrast to the uni- 
versal scaling of 0(iV 2//3 ) for most random networks 
with finite size percolation threshold. Interestingly, 
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for r approaching 3, where the network assumes a fi- 
nite size percolation threshold, the universal scaling 
exponent 2/3 is recovered. On the other hand, for 
p < 1, the size of the largest connected component 
scales at most as fast as JV( T-2 '/ T . For fixed but large 
N, we find that the size of the largest component as 
a function of p for all p, and we recover the results 
in for 1 -c p -C A 1 / 7 ". In general, for large p 1 
and N, the size of the largest connected component 
is shown to scale as ~ p 1 / i^- T > J<f 1 ~ 1 / T . 



C. Applications 

The results in this paper also have applications to 
peer-to-peer (P2P) data mining or search algorithms 
on unstructured complex networks. Percolation at a 
multiple of the threshold is at the heart of a scal- 
able P2P search algorithm, called percolation search 
algorithm, introduced by the authors in Percola- 
tion search uses a probabilistic broadcast algorithm 
to reduce the number of communications necessary 
for finding contents in a large scale complex network 
that has a heavy-tailed power-law degree distribu- 
tion. The percolation probability q will correspond 
to the probability with which nodes of the network 
communicate a message to any of their neighbors. For 
the search algorithm to work efficiently, one needs to 
choose a value of q that is as small as possible and 
yet results in a large enough connected component; 
that is, choose q to be a constant multiple of the 
threshold q c . The scaling of this connected compo- 
nent directly specifies the scaling of the P2P search 
traffic with network size. The results in this paper 
form the theoretical basis for showing that this traf- 
fic scales sub-linearly with the network size for any 
heavy-tailed random power-law network. 

The rest of this paper is organized as follows: In 
Section [n] we briefly review the generating function 
formalism |lj for calculating the percolation prop- 
erties of random networks on given degree distribu- 
tions. We then specialize the approach to the case of 
percolation at a multiple of the threshold for heavy- 
tailed random power-law networks and derive the 
scaling relation of the size of the connected com- 
ponent as a function of both q/q c and the network 
size. A number of simulations are reported in Sec- 
tioii lIIII which verify our results. Concluding remarks 
are made in Section ITVl 



II. CALCULATION OF THE CLUSTER 
SIZES 

In this section, we introduce the analytic approach 
of the paper. 



A. Modeling Finite Size Effects 

Consider a random graph with degree distribution 
Pk, k = 1, 2, k max . Thus, the probability that a 
randomly chosen node of this graph has degree k is 
Pk, where 1 < k < k max . We assume that the max- 
imum degree is such that every network is expected 
to have at least one node of that degree: Npk max > 1. 
For a power-law random graph pk — A T k~ T for 

h max 

), we have: 

A T « CO") (1) 

kmax ~ ■, / (2) 
Ar' 

This choice of k max reflects the maximum scaling of 
kmax that results in the expected number of high de- 
gree nodes to be only one. 

Thus, as N increases , k m ax will increase and the 
variance of the distribution will become unbounded 
for t < 3. We focus on the parameter regime, 
2 < t < 3, where the networks have degree distri- 
butions with unbounded variance but bounded mean 
(except at r = 2, where it increases only logarith- 
mically). Though we are interested in finite size ef- 
fects, we focus on the case of large values of N where 
(k 2 ) » (k). 

An accompanying assumption we make is that the 
generating function approach, which is meant to deal 
with infinite size networks, is still applicable to ran- 
dom power-law networks of finite, but large size, N. 
The same assumption is implicitly made in |4j , which 
also uses the same generating function formalism. A 
formal treatment of the validity of this assumption 
is beyond the scope of this paper. It is not, how- 
ever, difficult to see that the scaling of the variance 
(and, hence, of k max ) with N will play an impor- 
tant role in any such formal approach, and for a dis- 
cussion on choosing k ma x ~ A 1 / 1 ", see 0, 0- We 
have verified the predictions made by the generating 
function approach via large-scale simulation results, 
some of which are presented in Section ITTT1 and in Ap- 
pendix^] The simulation results in AppcndixlFllsliow 
that the predicted infinite-size percolation thresholds 
are matched closely by numerical estimates obtained 
from synthetic finite-size networks; the match is very 
good even for relatively small size networks. 

B. Percolation as a Branching Process 

Just above the criticality, the largest component of 
the network is a tree. Percolation on a tree can be 
viewed as a branching process. Using this fact, we 
will treat percolation as branching process using the 
generating functions method introduced in 0. The 
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generating function for the degree of a randomly cho- 
sen node can be defined as Gq(x) — J2k=T Pk x>: - The 
generating function for the degree of a node arrived 
at by following one end of a randomly chosen link is 
G 1 (x) = G' (x)/G' (l). 

The bond percolation with probability q on a given 
graph is as follows: For each link of the network, 
delete the link with probability 1 — q and retain it 
with probability q, independently. The distribution 
of the size of the connected components of a random 
graph after bond percolation with probability q can 
in principle be found. Throughout the rest of the 
paper, we use results in 0,0]. 

Let u be the probability that a random link does 
not lead to an infinite set of nodes, one then has: 



1 



qGt 



(3) 



This equation is understood as follows: Take any 
random link. The probability that the link is deleted 
(and thus does not lead to a giant component) is 1 — q 
(hence the first term on the right hand side). Now 
follow the link to one of its random ends to arrive 
at a node V. G±(x) is the generating function of 
the distribution of the number of links that go out of 
V. Being a random graph, the probability of each of 
those links to lead into a giant connected component 
is again u. The probability that none of the remain- 
ing links of V go to an infinite component is G\(u). 
Hence we have (Eqn. |3J). 

As far as the nodes are concerned, the probability 
that a random node does not lead to an infinite num- 
ber of nodes is thus Go(u), and hence the fraction of 
nodes in the infinitely large component is: 



5 = 1- G (u) 



(4) 



Below the percolation threshold, the only solution 
to (Eqn. |3J| is u = 1. In fact, this can be used to 
show that the percolation threshold of any random 
graph on a given degree distribution is given by: 



= G' 1 (1)- 1 = 



(k) 



(k 2 ) ~ (k) 



(5) 



When working just above the percolation threshold 
(as is the subject of this paper), one expects that 
u < 1, and hence we look for solutions of the form 
u = 1 — 5 with 5 « 0: 



■5 
5 
(1 



l-q + qG x {u) 
l-q + qG x {l-8) 

l-Gi(l-tf) 



(6) 



Follow a particular edge, define S the probability 
that it leads to an infinite component. In order to 
solve for 6, one can write a Taylor series expansion of 
Gi(l-S): 



OO 



G 1 (l-6) = }^G{ n \l)- 

71=0 



X-sy 



(7) 



For heavy tailed random graphs, all these moments 
(except for the mean G (l)) are large (in fact infinite 
as N goes to infinity). Nevertheless, for any finite N, 
these moments are still finite and the expansions can 
be carried out. In this work, we deal with heavy-tail 
power-law random graphs. For these graphs we have 

(fc" +1 > « (»rr_?) k max {k n ). We use this property 
to make the following approximation: 



G\ n \l) 



(k) 



(k 



n+l\ 



(k) 



Finally, by approximating the sum as an integral 
we can get that for n > 1 



(k 



n+l\ 



A T h, 



n+2-r 

max 



(8) 



While these approximations may seem crude, we 
show in Section II I II that the predictions we make 
match the simulation results. 



C. Solving for the Scaling Exponent 
We can put (Eqn. |SJ| into (Eqn. 0): 

Gx(l-J) = f) G W(l)t^ 
(k n+1 ) (S) n 



n=0 
oo 

E 

n=0 



(k) 



+ ^ (k) (n + 2-r) n! 

A h 2 ~ T J2. (— Ah \ n 



(k) £j„!(„ + 2-t) 

Defining a finite size percolation threshold like any 
other random network as q c = l/G[(l), one gets: 



A k 2 ~ T 



(k) 



(3-r)(fc 2 ) 
3-r 

fcmaxQc 
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which gives: 



Or equivalently using (Eqn. I14fl : 



Gxil-S) « l + i— if] \, 6k T r , 
k max q c ^ x n\{n + 2-T) 

Inserting back into (Eqn. EJ: 



3-t 



q(3 - t) 



(Skr, 



kmaxQc £r{ n\(n + 2 - r) 

V unj maxJ 



;/' - 2 - r) 

n— 1 x 7 



For each constant r, define a 



and 2 



q(3-r) 

Sk max to get an equation that is independent of k max 
in the limit of large k max : 



az 



^ nil n + 2 - r) 

n— 1 v ' 



0) 



(Eqn. |nj is our first result of this paper. It clearly 
states that as far as fc mal is large, z would be inde- 
pendent of k max and hence N. Thus, the fraction of 
links in the giant connected component will scale as 
S = zk^ ax for some constant z depending only on r 
and q/q c . 

The number of links in the largest connected com- 
ponent after percolation, will therefore scale as: 

g(r,p) i (k)6N - {k)k m \ x zN (10) 

where z depends only on r, p and k max depends only 
on N, r. The above relation is valid as far as the ratio 

p = q/q c for 

qc « ( k )/(k*) « (k)k;- a l = (fc)A 1 -^ (11) 

When 2 < r < 3, (Eqn. [^J) can be simplified in 
terms of incomplete gamma functions. The incom- 
plete gamma function is defined for 5R(a) > as: 



T(a,z) 



x a - L e~ x dx 



(12) 



At z = it takes on the familiar value: T(a, 0) = 
r(a). For all z, the following holds true: 

T(a,z) = ~(T(a+ l,z) - e~ z z a ) (13) 
a 

The series expansion for T(a, z) is as below: 

r(a,z) = r(a)-z a J2 u Z i ^ ( 14 ) 

A — ' n\[n 4- a) 



~ (-z)» _ r( a )-r( a ,z) l 

^ n\(n + a) z a a 

Applying this identity to (Eqn. ^ with a = 2 — r: 
r(2 - r) - r(2 - r, z) 1 



2-T 



(16) 



Applying the above to (Eqn. [SJ, with a = 2 — r, and 
setting /? = ^§Ef 

/ 9g = r ( 3 - T ^3-r > ,) +e _,_ 1 

when 2 < t < 3. 

To see how the scaling results of 4] follow from 
117|l . lets consider the case z ^> 1. In this case, 
r(a, z) —> 0, hence: 

1 



-az w z T " 2 r(2-T) 

v ; 2 - r 

/3z « (2-T)z T " 2 r(2-r) - 1 

= z r - 2 r(3-r)-l 

If z T_2 r(3 — t) 3> 1, we can neglect 1 to obtain: 

0z w z T - 2 r(3-r) 

T(3-r)^ 1/(3 " r) 



g(3 - r)r(3 - r) 



l/(3-T) 



9c(r - 2) 

/rr4 ^\V(3-r) 

1/(3-t) / 1 I 4 r J 



V (t-2) 



In other words, the fraction of the links in the 
largest connected component 5 scales as: 



when 1 <C p <C k max and 2 < t < 3 and when the 
scaling of k max is chosen according to Eqn. [21 (i.e., 
A 1 / 7 "). The number of links in the largest connected 
component will therefore scale as 



(18) 



Similar result follows for r = 2 (see Appendix 
IA"ll . In particular, the fraction of links in the largest 
connected component can be shown to scale as ~ 
N~ x / 2 pin p and therefore, the total number of links 
in the largest connected component scales as 



(19) 



G(p, 2) - plnpVNlnN 



G 



Another interesting regime is where p — ► 1 + or 
z -C 1. This will correspond to percolation very close 
to the threshold q c . If z -c 1 we can look at the first 
few terms of (Eqn. 

(-*)" 



^ rc!(ra + 2 - r) 

n— 1 v 7 



2 



3-r 

-l 



2(4 



4- r i 

z^max = 2 (1 — a(3 — t)) - -k ma . 



2 1 



4-T 

3-r" 



= 2 |— ^((p-l)//^" 1 /* 
o — r 



In other words, the scaling with the order parame- 
ter is again linear close to the finite size percolation 
threshold. 



D. Links in the Largest Component After 
Percolation 



The parameter 8 is the fraction of total links before 
percolation that are in the largest connected compo- 
nent after percolation. With percolation threshold 
being small, most of the links are deleted during the 
percolation process. We are interested in the proba- 
bility that a link not deleted during the percolation 
process is part of the largest connected component. 
Equivalently, what fraction of the links that remain 
after percolation are part of the largest connected 
component? 

To answer this, note that the number of links that 
remain after the percolation is closely approximated 
by 

H(N,p)=pq c (k) 

Using (Eqn. Ill|) to calculate q c and the scaling re- 
sults for 6, the probability that a link left after the 
percolation is part of the largest connected compo- 
nent can be calculated as follows. 



G(N,p) 
H(N,p) 



plnp t = 2 

p l/(3-T) N 2/r-2 2 < T < 3 



(20) 



Interestingly, for r > 2, the fraction of links that 
are part of the giant connected component constitute 
to an infinitesimally small fraction of all links after 
the percolation. Only for r = 2 a finite fraction of all 
the links that remain after the percolation belong to 
a single connected component. 



III. SIMULATIONS 

Figure [21 shows the scaling of the largest compo- 
nent after bond percolation as a function of the size 
of the original network N when q/q c is a constant @ 
using network sizes ranging from 1, 000 — 1, 000, 000 
nodes for different values of r. We compare the sim- 
ulation result to the scaling we would predict (i.e., 
TV 1- "). Simulation results suggest that our scaling 
law is correct as long as k ma x is very large compared 
to z and q/q c is in the order of one. 



IV. CONCLUDING REMARKS 

We investigated the properties of the size of the gi- 
ant connected component just above the percolation 
threshold in heavy-tailed power-law random graphs, 
for which the percolation probability is known to be 
vanishingly small. By normalizing the percolation 
probability by the percolation threshold, we were able 
to trace the scaling behavior of the size of the giant 
connected component at very small percolation prob- 
abilities. In particular, we showed that <5, the frac- 
tion of links in the giant connected component close 
to the percolation threshold, is proportional to the 



factor 



for kn 



A. High Degree Nodes and the Giant 
Connected Component 

Let us address the question of which nodes are most 
likely to be in the giant connected component when 
,r-2 

p = p is in the order of one. Note that the 

3 — T 

percolation probability is q c = (k)/((k 2 ) — (k)) oc 
k^ax for large k max , when 2 < r < 3. As such, the 
probability that any node with k links will have any 
edges left after percolation (with f) in the order of 
one) is around ~ kk^^,. For this probability not to 
be negligible, we must have k ~ k^J x or greater. Put 
in other words, take any node with degree kf r ^ a r x ~ e for 
any finite e > 0, then with probability one this node 
will lose all its edges after the bond percolation and 
will not be in the largest connected component after 
percolation at a multiple of the threshold. 



B. Percolation Search 

The percolation search algorithm, developed by 
the authors, is based on a probabilistic broadcast- 
ing scheme, with as small a probability of broadcast 
as possible. The success of this algorithm depends on 
finding a percolation probability which would result 
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q = 1 .1 q_c, tau=2.0 
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(our prediction) 
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q = 2.0 q_c, tau=2.5 
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10000 
nodes 



FIG. 2: Examples of the scaling of the size of the largest connected component with the network size N for different 
values of r and q/q c - From top down, r = 2.0, q/q c = 1.1, r = 2.2, q/q c = 20.0 and r = 2.5, q/q c = 10.0. The scaling 
AT 1 - 1 / 1 " is depicted on all the plots. Even in the case of the top plot, where the largest component is as small as 100 
our scaling predictions provide good matches to the simulations. 
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in most of the high-degree nodes to fall into one giant 
connected component, while as few low- connectivity 
nodes would be present in the same component. This 
will ensure that most of the search traffic will be car- 
ried out by high capability nodes that have assumed 
large degrees, while the low connectivity nodes will 
only occasionally participate in a search. Moreover, 



since a query message will be passed only along a few 
edges, the protocol will result in low overall traffic. 
The results in this paper (see the preceding discus- 
sions) show that random power-law networks are ide- 
ally suited for percolation search, and broadcasting 
with probability just above the percolation threshold 
leads to high query hit rates. For more details see p|. 
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APPENDIX A: CASE OF t — 2 

In the case of r = 2, (Eqn. reads as: 
\- {-z) n 



percolation threshold of finite networks tau=2 



nin 



(Al) 



The above series is related to the Ei(z), or expo- 
nential integral function: 

/oc — t 
— alt . 
-z * 

The series representation of Ei(z) is: 

°° z n 
run 

n=l 

Where 7 is the Euler-Mascheroni constant and is ap- 
proximately 0.5772. This gives us: 

00 / \ n 

J2 [ -^-=Ei(-z)-j-hx\z\ (A2) 

n—l 



We can then write (Eqn. I Al|) as: 

— az — Ei(—z) — 7 — In z 



(A3) 



when r = 2. For z ^> 1, Ei(—z) — ► 0, or z/hx{z) 
cv _1 = p or equivalently, z « phi(p). 
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FIG. 3: For a finite network we define S(q) as the size of 
the largest component as function of percolation proba- 
bility q and find the value of q where dS/dq is maximized 
(marked by the Simulation points). We compare that 
value with the percolation threshold of an infinite network 
with the same degree distribution: q c — {k)/({k 2 ) — (k)). 
We see that the ratio of these two values is approximately 
constant and tends towards unity as the network grows 
in size. 



APPENDIX B: VERIFYING THE 
PERCOLATION THRESHOLD 



The discussions in this paper highlighted the sig- 
nificance of q c (N) = (k)/{{k 2 ) — (k)) as the true 
finite size percolation threshold. Strictly speaking, 
however, q c {N) > (for any fixed N) is the per- 
colation threshold of an infinite size random graph 
whose degree distribution is given by the truncated 
PL: p k = A T k~ T for k G (l,k max ), and p k = for 
k > k max (we let k max scale as N 1 ^ , and hence 
is fixed for a given N). We know from percolation 
theory 0, Q that for a random graph with a fixed 
kmax, and the degree distribution satisfying certain 
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additional technical conditions, the infinite size per- 
colation threshold can be shown to be the limit of a 
uniformly convergent series, and that the phase tran- 
sition for any network of finite but large size, M, 
will also happen close to the infinite size percolation 
threshold. The exact dependence of the network size 
Mq (above which the infinite size percolation thresh- 
old will be a good approximation for finite size perco- 
lation) on the nature of the degree distribution and 
kmax (especially when it is large) is currently un- 
known. In this paper, we are assuming that for any 
fixed but large N (which determines k max , and hence 



q c (N), for the family of random networks), the infi- 
nite size percolation threshold is also a good approx- 
imation for a finite network of size N. 

Figure 13 directly compares what we might call the 
empirical threshold to the infinite network model. If 
S(q) is the size of the largest component for a fixed 
network size with a bond percolation probability q, 
wc define the empirical threshold as the point where 
dS(q)/dq is maximized. For all networks, the em- 
pirical threshold is very close to the infinite network 
model, and as N grows, these two agree more closely. 



